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Abstract. We use models of magnetized hydrogen atmo- 
spheres to fit the soft X-ray thermal spectrum of the Vela 
pulsar obtained by ROSAT. The distance and hydrogen 
column density required by our fits are in good agreement 
with other independents estimates and can be considered 
as an argument for the presence of hydrogen at the surface 
of this pulsar. 

The low temperature obtained, T e °° = 7.85 ± 0.25 x 
10 5 K, is below the predictions of the 'standard' model of 
neutron star cooling and strengthens previous claims for 
the necessity of the occurrence of fast neutrino cooling in 
this neutron star. 



1. INTRODUCTION 

The Vela pulsar is the youngest pulsar for which there is 
strong evidence that thermal emission from the surface 
of the neutron star has been detected (Ogelman, Finley 
& Zimmermann 1993; Ogelman 1995). A reliable mea- 
surement of its surface temperature based on modeling of 
the surface thermal emission is of outmost importance for 
comparison with models of neutron star cooling and could 
give us evidence for the past occurrence of fast neutrino 
emission. Moreover, since the chemical composition of a 
pulsar's surface and its state (gaseous, liquid or solid) are 
still unknown, the 'successful' fit of the observations can 
also give us invaluable information about this question. 

2. MAGNETIZED HYDROGEN ATMOSPHERE 

There are many possibilities for the structure of a pul- 
sar's surface and a magnetized hydrogen atmosphere is 
the most 'natural' and simplest choice although it has a 
priori no strict reason to be accepted (e.g., where does the 
hydrogen come from ?). Assuming the presence of an at- 
mosphere, the magnetic field has three main effects: 1) the 
opacity depends strongly on the magnetic field direction 
and on the photon polarization, 2) for the polarization giv- 
ing the main contribution to the emergent flux it is much 
smaller than in the nonmagnetic case and its dependence 
on the photon energy is weaker, 3) the ionization energy 
is strongly increased and at low T e an absorption edge ap- 
pears within the range of the ROSAT PSPC. The problem 
of ionization equilibrium in a magnetized hydrogen atmo- 
sphere is not yet completely solved but at T e > 10 6 K the 
whole atmosphere is fully ionized, so its structure can be 
calculated accurately and the spectra we are using are re- 
liable (Shibanov et al. 1992; Pavlov est al. 1994, 1995). 



3. SPECTRAL FITS 

We performed a % 2 search of the whole parameter space 
using the ROSAT data as presented by Ogelman, Finley 
& Zimmermann (1993), restricting ourselves to E < 1.5 
keV to eliminate the hard tail. Due to the high tempera- 
ture, several simplifications occur when considering emis- 
sion within the ROSAT PSPC band: 1) the emitted flux is 
weakly dependent on the field strength at B > 10 12 G, 2) 
the detectable flux practically depends only on the mean 
'effective' surface temperature, not on the actual distri- 
bution of the temperature along the surface (Page 1995a, 
Page & Sarmiento 1995), 3) the red-shift is equivalent to 
a change (red-shift) of temperature and 4) the model de- 
pends only weakly on the surface gravity. 

As a first approximation we are thus able to do spec- 
tral fits with only three parameters: effective temperature 
at infinity, T e °°, distance D and column density Njj. We 
assume a mass of 1.4 M and a radius of R = 10 km, 
i.e., R°° — 13.6 km, but since the received flux scales 
like (R°° /D) 2 our results can be extrapolated to other 
radii with reasonable confidence. The final result of our 
X 2 search is shown in Figure 1. 

The pulsations of the soft X-ray emission from Vela 
show that the temperature of the polar cap region may 
be even larger than the mean effective surface tempeature 
found in our fit (T^/ ~ 9 x 10 5 K ) and the fully ion- 
ized atmosphere fit might be quite reliable for this object. 
However, at T^f < 10 6 K an absorption edge (not taken 
into account in the spectra we used) should appear: as a 
result our fits give too low a T e and our lower range would 
probably be pushed to higher values. 

4. DISCUSSION 

We obtain in our best fit a x 2 per degree of freedom of 
1.074; due to the low statistics any \ 2 P-d.f. of order one 
is acceptable. So the relevance of the model can only be 
assessed by discussion of the resulting values of the pa- 
rameters. 

Column density: N H = 2.35 ± 0.8 x id 20 crrT 2 . 
Other components of the Vela PSR + SNR give in- 
dependent estimates of Njj- Modeling of the pulsar 
'hard tail' soft X-ray emission needs values of 4.3 ± 1.0 
or 4.4 ± 1.5 xlO 20 cm" 2 and modeling of the sur- 
rounding compact nebula gives values of 2.0 ± 0.5 or 
3.0 ± 0.5 xlO 20 cm~ 2 (Ogelman et al. 1993) while the six 
'bullets' detected in the SNR give N H 's between 2.4 ± 1.9 
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Table 1. Reference objects for interstellar absorption 
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Fig. 1. Confidence contours for the three parameters (T e , Nh 
and D), projected onto the T e — Nn plane. The grey area shows 
the 'lcr' lower limit on the distance from the pulsar dispersion 
measure. The dot shows our best fit: T e °° = 7.04 x 10 5 K, 
(Tsurf = 9.2 x 10 5 K), Nh = 2.35 x 10 20 cm -2 and D = 285 pc. 



and 5.4 ± 2.1 x 10 20 cm" 2 (Aschenbach et al. 1995). In 
addition, Table 1 list several reference objects in the line 
of sight of the Vela pulsar for comparison with our re- 
sults which show that Njj reaches 10 21 cm -2 well above 
500 pc and probably needs at least 100 - 200 pc to reach 
10 20 cm -2 . In short, a Nh between 1.5 and 3 x 10 20 cm -2 
seems reasonable. 

Distance: D = 300 ±120 pc. The pulsar's distance of 500 
± 125 pc (Taylor et al. 1993), confirmed by interstellar 
scintillation measurements (Gupta 1995), is compatible, 
in its lower range, with our deduced value, in its upper 
range. Moreover, increasing the star's radius R°° increases 
the fitted distance D in the same amount. Our results may 
possibly favor a larger radius: R = 12 km (i?°° = 15.5 km) 
would increase D by 14%. 

Our 68% confidence range ('ler') combined with the 
pulsar's distance (at 'ler', the grey area in Fig. 1) thus 
implies that T e °° = 7.85 ±0.25 x 10 5 K (for a 10 km radius) 
and N H = 1.80 ± 0.25 x 10 20 cm" 2 . The corresponding 
surface temperature T^f is above 10 6 K where our atmo- 
sphere models are most reliable. 

Hydrogen atmosphere ? The consistency of our result for 
D and Nh with values obtained with other methods is 
an argument in favor of the presence of a hot (~ 10 6 K) 
magnetized ( > 10 12 G) hydrogen plasma at the surface of 
the Vela pulsar in sufficient amount to provide an optical 
depth of unity, i.e., at least a few grams per cm 2 . 
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20.48 ± 0.10 
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-3.6 


20.84 ± 0.10 
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-3.1 


21.20 ± 0.05 
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a: Fruscionc. A. et al. 1994. 
b: This work 

c: Ogelman, Finlcy. Zimmcrmann 1993. T° 
fit. 



1.6 X 10° K, blackbody 



5. COOLING 

The Vela pulsar has already been claimed several times to 
be cooler than the predictions of the 'standard' model of 
neutron star cooling and our results reinforce such claims. 
If we accept our deduced T e and believe neutron star cool- 
ing models then our results require fast neutrino emission 
by, e.g., direct Urea processes, kaon or pion condensates, 
which has been supressed by baryon pairing (e.g., Page & 
Applegate 1992). A critical temperature T c of the order of 
1 — 3 x 10 9 K for pairing in the inner core is then neces- 
sary, independently of the exact fast cooling agent (Page 
1995b). 
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